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Abstract 

The Higgs boson H of the Standard Model can be searched for in the channels pp/pp —> 
ttH + X at the Tevatron and the LHC. The cross sections for these processes and the 
final-state distributions of the Higgs boson and top quarks are presented at next-to- 
leading order QCD. To calculate these QCD corrections, a special calculational technique 
for pentagon diagrams has been developed and the dipole subtraction formalism has been 
adopted for massive particles. The impact of the corrections on the total cross sections is 
characterized by K factors, the ratios of the cross sections in next-to-leading order over 
leading order QCD. At the central scale po = (2m t + Mjj)/ 2 the K factors are found to 
be slightly below unity for the Tevatron (K ~ 0.8) and slightly above unity for the LHC 
(K 1.2). Including the corrections significantly stabilizes the theoretical predictions for 
total cross sections and for the distributions in rapidity and transverse momentum of the 
Higgs boson and top quarks. 
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1 Introduction 


Even though the Standard Model (SM) has been tremendously successful in describ¬ 
ing matter and forces in particle physics, the third component of the model, the Higgs 
mechanism [ [3], has not been established yet experimentally. This component, however, 
is crucial for the closure of the SM in a mathematically consistent form, thereby allowing 
for systematic perturbative expansions and precise predictions in the electroweak sector [ 
. The search for Higgs bosons [ |3j. [|] is therefore one of the most important experimental 
programs in present-day high-energy physics [ || . 

For a light Higgs boson the electroweak interactions can be extended perturbatively up 
to the grand unification scale. The validity of this extrapolation is backed by the fact that 
it provides a qualitatively correct estimate of the electroweak mixing angle. Demanding, 
therefore, the perturbative extrapolation to hold up to the grand unification scale, the 
mass of the Higgs boson is bounded to be below 180 GeV in the SM [ An upper limit 
in the same range is also obtained experimentally by evaluating radiative corrections to 
electroweak precision observables, i.e. Mh < 196 GeV at the 95% confidence level [[?J. A 
lower experimental mass limit of 114.4 GeV has been set by the LEP experiments [ ||. 

In the near future, the search for Higgs bosons will continue at hadron colliders, the 
proton-antiproton collider Tevatron [ 3j with a centre-of-mass energy of 2 TeV, followed 
by the proton-proton collider LHC [ []J with a centre-of-mass energy of 14 TeV. Various 
channels can be exploited at hadron colliders to search for Higgs bosons in the intermediate 
mass range. Among these channels, Higgs-boson radiation off heavy top quarks [ |9j. [TI]. |TT[ 
plays an important role: 

pp/pp —> ttH + X. (1.1) 

Although the expected rate is low at the Tevatron, a sample of a few but very clean 
events could be observed for Higgs masses below 140 GeV [ |T2| . At the LHC, Higgs- 
boson radiation off the top quarks is an important search channel for Higgs masses below 
~ 125 GeV [ [HI [bl. Moreover, analyzing the ttH production rate at the LHC can pro¬ 


vide valuable information on the top-Higgs Yukawa coupling in relation to other Yukawa 
couplings [ [15], [L6|] . Precision measurements of the absolute value of the ttH coupling can 
be completed later at e + e _ colliders [ [L7|, [18], [19], [20| [21], [22[ . 


Theoretical predictions for cross sections that are based merely on leading-order (LO) 
QCD are notoriously imprecise. They are plagued by considerable uncertainties owing 
to the strong dependence on the renormalization and factorization scales, introduced by 
the QCD coupling and the parton densities. Therefore, higher-order QCD corrections are 
needed for satisfactory theoretical predictions. 

The cross sections for the processes (Q are found in next-to-leading order (NLO) 
QCD by convoluting the cross sections da ab of the parton subprocesses ab —> ttH ( +k ) 
with the parton distributions f a and fb of the initial-state hadrons, both evaluated in 
next-to-leading order: 


da(hih 2 -*■ ttH + X) = J2 f d £« /o'te*; /4) J d 6> fb' 2 (fb] /4) dcr ab (s = f a f 6 s; /4, /4)- 

a,b=q,q,g 

( 1 - 2 ) 

In this convolution the partons a and b carry fractions £ a and of the original momenta 
of the hadrons h\ ( =p) and h 2 ( —p/p ), respectively. The total hadronic centre-of-mass 
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energy is denoted by yfs, whereas the partonic subprocess energy is given by V§. The 
factorization and renormalization scales are denoted by p,p and /i#, respectively. In general 
these two scales will be identified. 

A brief report on the radiative QCD corrections for the processes ( |1.1|) in NLO QCD 
has been presented for both the Tevatron and the LHC in Ref. [0- The results of that 
study were found to be in agreement with a parallel investigation for the Tevatron in 
Ref. [0], which was based exclusively on the dominant partonic process qq —> ttH. As 
expected, in complete calculations of the NLO QCD corrections the scale dependence 
is reduced significantly, and stable theoretical predictions can be derived for the cross 
sections. 

At the technical level there are two main obstacles in calculating the NLO QCD cor¬ 
rections. First, the one-loop pentagon (5-point) diagrams, involving both soft/collinear 
singularities and massive particles, have to be calculated in dimensional regularization. 
To this end, we have developed a calculational technique, exploiting the fact that the 
singularity structure of the pentagon diagrams can be derived in a universal way. This 
enables us to switch to regularization schemes that are more suitable for analyzing pen¬ 
tagon diagrams, without losing regularization-scheme-dependent constant terms. Second, 
the extraction of the singularities in the real part of the NLO QCD corrections has to 
be performed in a numerically stable way. For this purpose we have adopted the dipole 

[This is the first application of 


subtraction formalism [ ^3(] for massive particles [ ^4 


this particular method to a complex NLO QCD calculation with massive quarks.] As an 
additional cross-check, the partial results for qq —> ttH have been recalculated by means 
of the phase-space slicing method. 

Besides the total cross sections, final-state distributions of the Higgs particle and the 
top quarks in rapidity y and transverse momentum pt will be presented. As in the 
case of the total cross sections, including the NLO corrections significantly reduces the 
scale dependence of these distributions. It should be noted that based on this complete 
calculation, the NLO QCD corrections to arbitrary distributions can be investigated. 

The report is divided into four sections. In the next section we describe the calculation 
of the NLO QCD corrections at the partonic level. In the third section the hadronic 
results for the SM reactions pp/pp —> ttH are presented, comprising the total hadronic 
cross sections as well as the distributions in rapidity and transverse momentum. The 
analytical and numerical results of our study are summarized in the last section. Finally, 
the appendices provide some useful formulas, such as scalar one-loop integrals, results on 
colour algebra, and splitting functions. 


2 Calculation of the NLO corrections 

2.1 LO cross sections and conventions 

The LO hadronic processes pp/pp —> ttH + X proceed at the parton level via qq 
annihilation and gg fusion: 

?(Pi) + ?(P2) t(p 3 ) + t(p 4 ) + H(p 5 ), (2.1) 

ff(Pi) + g(P2) -»■ t{p 3 ) + t(p 4 ) + H(p 5 ), (2.2) 
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where the momenta of the particles are given in brackets. The momenta obey the on-shell 
conditions p\ = p\ = 0, p\ = p\ = mf, and pi = Mjj. For later use, the following set of 
kinematical invariants is defined: 

s = ( Pl +p 2 ) 2 , 

Sij = {jPi +Pj) 2 , i,j = 3,4,5, 

Uj = (Pi-Pj) 2 , * — 1)2, j = 3,4,5. (2.3) 

A generic set of LO diagrams is shown in Fig. |1| all LO diagrams not shown differ only 
in the point where the Higgs line is attached to the top-quark line. 

We express the amplitudes J\/[ 99 ’ 99 of the LO and one-loop diagrams for the qq an¬ 
nihilation and the gg fusion channel in terms of colour factors C 99 ’ 99 , standard matrix 
elements (SME) A4 99 ’ 99 , and invariant functions F??’ 99 . The SME include the Dirac 
structure, spinors, and polarization vectors in a standard form, while the F 99 ’ 99 comprise 
all remaining factors, such as propagators, couplings, and standard one-loop integrals. 

To introduce a compact notation for the SME, the tensors 


r {?„...} = v q(P2) I 1 ) 7m, 7m7^, • • •} Uqipi), 

= MPa) {1, 7m, 7m7v, • • •} W(p 4 ), (2.4) 

are defined with obvious notations for the Dirac spinors u q (pi), etc. Furthermore, symbols 
like T p are used as shorthand for the contraction T^p 9 . For the qq channel we introduce 
the following set of SME, 


KA qq — F?9’M rtt 

J 1 { 1 ,2,3,4} L L {m,MPi>MP2,MPiP2}’ 

A A 99 — FW.MP3P4 pit 

J l {5,6,7,8} 1 L (MiMPiiMP2,MPiP2)’ 

M 99 _ pqq,P3 rtt 

J 1 {9,10,11,12} 1 1 {1,P1,P2,P1P2}’ 

a a 99 — rqq-.P'i rtt 

J l {13,14,15,16} 1 1 {l,Pi,P2,PiP2}’ 

a A 99 — rqq,P"P3 rtt 

J 1 {17,18,19,20} {m^m^Pi :M l/ P2>M t 'PiP2} ’ 

Ayf 99 — rqq-jj^TM rtt 

J 1 {21,22,23,24} 1 1 {m^MMPi ,MMP2,MMPiP2} ’ 

a A q9 — rqq-i ii/ p rtt 

J 1 {25,26,27,28} {MJ'P,Mt'PPl,M I 'PP2,M I 'PPlP2r 


(2.5) 


The LO amplitude A4l 9 0 °f Hie QQ channel is proportional to the colour factor 


C 99 = \ c ® A c , 


( 2 . 6 ) 


where the hrst Gell-Mann matrix A c belongs to the qq state, the second matrix to the tt 
state. The one-loop amplitudes M 99 _ loop also involve other colour operators, such as 1 <E> 1 
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and (A a A fe ) (8) (A a A b ), but they can be reduced to multiples of C qq and tensors that are or¬ 
thogonal to C qq and therefore do not contribute in the interference Xboiours (^lo)*T- ioop• 
Thus, C qq is the only relevant colour structure at one loop, and the produced tt pair is in 
a pure colour-octet state in the qq channel. Both the LO and the (relevant part of the) 
one-loop amplitude can be written as 


28 


M qq = C qq A qq , A qq = Y Fi q MT , 


i=1 


leading in particular to 


£ |A<r 0 | 2 = 32E |-4g,| 

spins, colours spins 


Explicitly, the LO amplitude can be cast in the form 


(2.7) 


( 2 . 8 ) 


Mfi, = - C™ 


a Tf&sQttH 


(Mf + Mf + 2 Mil) 


+ 


s 35 - m t 

1 

s 45 - m\ 


(.Mf + Mf - 2 M 


(2.9) 


from which the invariant functions Q i can be read off easily. The factor g tt n = m t /v 
denotes the SM Yukawa coupling strength, where v is the vacuum-expectation value of 
the Higgs held, and a s = g^/( 47t) is the QCD coupling constant. 

To dehne the SME for the gg channel, the polarization vector sf is introduced for each 
gluon with momentum pf Apart from the transversality condition e^pj = 0, the gauge 
conditions 

£iP2 = £2Pi = 0 (2.10) 

are used to simplify the algebra in the calculation of the amplitudes. Within this gauge, 
the following set of SME is complete, 


M 


Mf = 


■ /Vl {2,3,4} 

Mf 

A/I g9 

JVV {6,7,8} 

Mf 

A A99 

yVl {12,13} 

4,15} 

Mil 

7 ,18| 


99 

{19,20,21,22,23} 


pit _ pit _ pit 1 pit 

1 PlP2£l£2 1 P2Pl£l£2 1 Plp2£2£l ' 1 P2pl£2£l 


= r f x {(eie 2 ), (£iPs)(£ 2P3), {eipf&pf} , 


pit _ pit 

1 £l£2 1 £2£l 


pit _ 

1 P1P2 

pit 

1 Pl£l£2 


pit 
1 P2P1 


x {(ei£ 2 ), (MPsX^Ps), fm)(£ 2 Pf} , 


~ntt 

1 pi£2^1 


~ptt _ ~ptt 

1 PlP 2 £l 1 P2Pl£l 




x {(£ 2 Ps), (£2^4)}, 

p l6l X {(£2P 3 ),(£2P 4 )}, 


-ptt 


P 2£l 


X {{£2P3),(£2Pf} , 


~ntt 
1 P1P2 


-ptt 
1 P2P1 


(£iPs)(£2P4), 


= r ei X {(£2P3),{£2Pa)}i 


T“ x {(eie 2 ), (eip 3 )(e 2 P3), (£iP3)(£2Pa), (m^X^s), (eiPdfap*)} , 
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Mil = r?(£ip 3 )(£2P4), 



M 99 


1 

Pl^P2,Sl^£2 


i = 1, - - -, 24. 


(2.H) 


In order to exhibit the Bose symmetry, the SME have been divided into the two subsets 
M 99 and M 99 which transform into each other by interchanging the incoming gluons. 
Obviously the first eight SME of the two sets are Bose symmetric or antisymmetric, 


M 99 = +M 99 , i = 1,2,3,4, 

M 99 = - M 99 , i = 5,6, 7,8. 


Alternatively we may write 


MH +i = Ml L = M" 


l 8+i 


8-hi 




i = 1,..., 16, 


( 2 . 12 ) 


(2.13) 


for the independent crossed SME M 99 with i > 8. 

For the gg channel it is convenient to introduce the three colour operators 


C 99 = 5 C1C2 1, Cf = i/ ClC2C A c , Cl 9 = d ClC2C X c , (2.14) 


where c n is the colour index of gluon n, and the matrices 1 and A c act on the ti state. 
The totally antisymmetric and symmetric SU(3) structure constants f abc and d abc are 
defined in the usual form. Parts in amplitudes that are proportional to different C 99 do 
not interfere with each other owing to the orthogonality relations: 

Tr [c 99 ^C 9 k 9 ^j = c 99 Sjk with c{ 9 = 24, c 9 2 9 = 48, cf 9 = ^. (2.15) 

Obviously C 99 corresponds to a colour-singlet state of the ti system, while C 99 and C 99 
describe the two different octet states. Both the LO and the one-loop amplitude can be 
written as 

3 40 

M" = J2 C f A f 1 A f = E F ijM 9 i 9 , (2.16) 

3 =1 *=1 

leading in particular to 

y iMfoi 2 = y^y Wo/. (2.17) 

spins, colours j= 1 spins 

Explicitly, the LO amplitude reads 

■C 99 + C 99 + C 99 ^ (^M <\irect + “d7f us i on ^ 

C{‘ - Cf + CJ a ) (M„««« d - i M tolon ) (2.18) 

Mf + (4 + «13 - mA (X? + Mr) 


with 

•i Alirect 


fS/IK 


\A" — 

JVl LO ~ 


7T QsQttH 


rriiMt 1 T — m? 1 
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- AMI 9 + Mf - 2M\l - 2M" 


12 


na s g ttH 


{ti3 ~m 2 t )(t 2 4 - rnf) 


Mf 9 + |(Mf + Ml 9 ) + Ml 9 + 2M 13 


99 


4Mf - 2M" 


12 


+ 


QttH 

(■S 35 - "h 2 )(^4 - rn%) 


Mi 9 +{^ + t 2 t~m 2 t \(M 9 2 9 + M 9 5 9 ) 


AMf + Mf + 2M?f + 2M?f 


l 13 


-^crossed 

^fusion 


= M, 


direct » A gg »~ A gg "> 


9ttH 


s(s 45 - m 2 ) 

2t[ QL s (JtfH 


+ 


5^35 - wit ) L 


(^13 — ^23 )Mf + 2 Mf 
(t 2 4 — ti 4 )Mf + 2 Mf 


— — Mi 


fusion 


tu~t 2i ,M° g ~M? 3 


(2.19) 


The terms in Eq. ( |2.18|) proportional to M direct correspond to the three “direct graphs”, 
where two parallel gluons couple to the top-quark line (see 3rd graph in Fig. [I]). The terms 
proportional to M crossed result from the direct graphs by crossing, i.e. by interchanging the 
two gluons (see 4th graph in Fig. |lj). The terms proportional to Mf usion correspond to the 
gluon-fusion graphs (see 2nd graph in Fig. |T[) which do not receive crossed counterparts. 

Some comments on the SME and their evaluation ought to be added. The M 99 ’ 99 
defined above lead to a unique representation of all LO and one-loop diagrams, as indicated 
in Eqs. (|2.7|) and ( j2.16|) , if only the Dirac equation, the transversality of the polarization 
vectors, and the gauge conditions ( [2. 10| ) are used. The number of SME could actually be 
reduced further by exploiting discrete symmetries and by using the four-dimensionality 
of space-time, inducing relations among the SME given above. However, eliminating 
some of the SME does not simplify the numerical evaluation significantly. The explicit 
calculation of the SME was carried out in two different ways. First all interference terms 
Sspins (M 99,99 )* Mf’ 99 were calculated by making use of the general polarization sums for 
Dirac spinors and polarization vectors. In a second independent approach we directly 
expressed the M 99 ’ 99 in terms of Weyl-van der Waerden spinor products [ ^3] for each 
hclicity configuration. 

The final step in the calculation of the partonic cross sections and distributions in¬ 
volves the 4-dimensional integration over the final-state three-particle phase space. By 
introducing an intermediate “ff-state” with virtuality S34, one obtains 


d < 0 = 


_ l,a.b 

2 s fcav 


(27 t ) 4 8 S34 

x 0 ([a/§ - M h } 2 - S34) 


K{s ™' m >' m ‘ ] e( y| _ 2mi - M„) 0( S3 4 - 4m?) x 

8 S34 8 s 


E |Mlq| 2 ds 34 dflj dcos(6^ M ), 


spins, colours 


with 


n(x,y,z) = \J'x 2 + y' 2 + z 2 — 2xy — 2xz — 2yz . 


( 2 . 20 ) 


( 2 . 21 ) 
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Figure 2: Representative diagrams of the gauge-invariant subsets (i)-(iii) of one-loop QCD 
diagrams, as described in the text. 


The labels ab indicate the specific initial state of the partonic process. The solid angle 
of the top quark is defined in the rest frame of the “it-state”, using the direction of 
flight of the “tt-state” in the initial-state centre-of-mass frame as z-axis. The angle 
denotes the polar angle of the Higgs boson in the initial-state centre-of-mass frame, using 
the direction of the initial-state beams as z-axis. The corresponding azimuthal angle is 
trivially integrated to 27T in view of rotational invariance about the initial-state beams. 
The factor k results from the average over the initial-state spins and colours: 

1/19 — J_ U99 — _ U99 = urn — J_ (n OO'i 

av 36’ av 256’ av av 96 ’ [ ' 


where k™ and are given for future use. The differential cross section including virtual 
NLO corrections can be obtained in a straightforward way from Eq. ( 2.2C|) by substituting 
the correct matrix element. 


2.2 Virtual corrections 

2.2.1 One-loop diagrams and calculational framework 

Among the one-loop QCD diagrams, three different gauge-invariant subsets can be 
distinguished, for both the qq and gg channels. Representative diagrams of these subsets 
are shown in Fig. |2j. Subset (i) comprises all closed-quark-loop graphs where the Higgs 
boson couples to the external top-quark line, (ii) includes all graphs where the Higgs 
boson couples to a closed quark loop, and (iii) is formed by all other diagrams, which 
include at least one gluon in the loop. Counting closed quark loops only once per quark 
flavour, we get 2, 2, and 25 graphs of the respective subsets in the qq channel; and 6, 
20, and 108 graphs in the gg channel. The counting is performed in the Feynman gauge, 
omitting tadpole-like graphs and external self-energy corrections consistently. Of course, 
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the subsets (i) and (iii) receive counterterm contributions from renormalization, which 
are considered in Section \ 2.2.2\ in more detail; the second subset is not renormalized in 
NLO, since Hgg and Hggg couplings do not exist in LO. 

In NLO the one-loop amplitudes contribute via interference with the corre¬ 


sponding LO amplitudes A41q". Using the colour decomposition rules of Section ETTL the 


square of the sum of the LO amplitude and its virtual NLO correction reduces to 

E |Vg, + (i-, 0 «p,l 2 = 32 E (Wtl 2 + 2Re{KWM?l ]oop }) , 


spins, colours 


/U l-^LO+(l-loop)| 2 
spins, colours 


spins 


= E<f E (Mlo/ + 2Re{(^f o ,,)Mf_ l0opJ .}) .(2.23) 

j =1 spins 


The evaluation of the virtual NLO correction was performed in two completely in¬ 
dependent calculations, leading, in particular, to two independent computer codes. In 
the first evaluation, the Feynman graphs have been generated by using FeynArts [ |26| . 
These amplitudes have subsequently been reduced in terms of SME with the help of 
Mathematica. Following standard techniques of one-loop calculations [ [2?]. |28], [h]. |3(|, |3lJ , 
the tensor integrals are algebraically reduced to scalar integrals. The treatment of the 
so-called “pentagon 1 ' diagrams, involving five propagators in the loop, requires particu¬ 
lar care and extensions of the standard methods given in Refs. [ [27| [F]. [2fj|. [F]. |3lJ . The 
treatment of the pentagons is described in Section [2.2.3| in detail, where the results on the 
scalar 5-point functions are listed explicitly. The results on the divergent scalar 3- and 
4-point functions are collected in Appendix |S|. Both ultraviolet (UV) and soft/collinear 
infrared (IR) divergences are treated in dimensional regularization. The renormalization, 
which removes the UV divergences, is discussed in Section [2 .2. 2| . The explicit structure of 
the IR. singularities occurring in the virtual corrections is given in Section |2.2.5| , and the 


cancellation against their counterparts in the real corrections is described in Section |273 . 
The algebraic output of the first calculation has been implemented in Fortran for nu¬ 
merical evaluation. The second calculation basically follows the same strategy, with a 
conceptual difference that is described in Section |2.2.3| \[v} . The amplitudes have been 
generated by hand and reduced with Form [ [52] , and the numerical evaluation is carried 
out in C++. 

For later use, we introduce the following definitions for the scalar and tensor integrals 
in analogy to Refs. [ E7], |28|, |29l, |3C 


(27r/r) 2e r D {1, q^, q^qu, ■ ■ •} 


iW,...)(?i,m ».m.) = ++~J dV doii 
C{o, ( 0 ,pp,...}(si,fctn 0 , mi, m 2 ) = —Er J d- ■} 
D{o,^u,...}(qi,q2,q3,rn 0 , m 1 ,m 2 ,rn 3 ) = J d D q- 

v ( \ [ a d iUtLA 

E{o^,...}{qi,q 2 ,q 3 ,Q 4 ,mo,m 1 ,m 2 ,rn 3 ,m 4 ) = — J d q —, 


m* J dodid,2 

( 2ng,) 2e 

R 2 J dod^dz 

(2vr/i) 2e [ D {1 ,q^q^q v ,...} 


17T- 


with d 0 = q 2 — ml + iO, dj — (q + qj) 2 — m 2 + iO, j = 1,..., 4. (2.24) 


































Here +iO is the infinitesimal imaginary part as required by causality, D — 4 — 2e is the 
space-time dimension, and /j is the arbitrary reference scale of dimensional regularization, 
which we may identify with the renormalization scale. The UV and IR divergences appear 
as poles in e; it is convenient to use the abbreviations 




r(i + a 

e 



A 2 (h) 


r(i + e) 



(2.25) 


for these divergences. [The top-quark mass has been used to define a characteristic energy 
scale.] If we distinguish between UV and IR divergences, we explicitly write A RV,IR (p). 


2.2.2 Renormalization 

The renormalization is performed in the MS scheme with the top-quark mass defined on 
shell. The top quark is decoupled from the running of the strong coupling a s (/i). Denoting 
the bare top-quark mass and the bare strong coupling as m t ° and respectively, these 
conditions fix the renormalization parameters in the transformations m ( ° = m t + 8rn t and 
« s ° = a s + 5a s as follows, 


Sm t 

m t 

doi s 

(y, s 


a. 
37T l 


3A rv (/i) +4 


cx s 

2n 


3 


11 

T 


) A™ (m,) + 


(2.26) 


where Nf = 5 is the number of light flavours. The divergence A RV (m t ) = T(1 + e)(47r) £ /e 
is independent of /j and originates from light-quark and gluon loops. The term propor¬ 
tional to Ai' v (p), on the other hand, originates from the top-quark loop in the gluon 
self-energy that is subtracted at zero-momentum transfer. In this scheme the running of 
the coupling a s (p) is generated solely by the finite contributions of the light-quark and 
gluon loops, while the top-quark contribution is absorbed completely in the renormaliza¬ 
tion condition and thus decouples effectively. The transformation for m t ° also fixes the 
renormalization of the Yukawa coupling, g^ tH = m^/v = (rn t + Sm t )/v, since v does not 
receive NLO QCD corrections. 

We renormalize the fields of the gluons, G a of the light quarks, ifjq, and of the top 
quark, all in the on-shell scheme, i.e. the wave-function renormalization constants 
5Zc, q ,t, defined by the transformations 


& a , M = (1 + SZ G / 2) G^, 


< = (1 + 8ZJ 2) $ = (1 + SZ t /2) A, (2.27) 


are adjusted to cancel the external self-energy corrections exactly. Distinguishing between 
divergences of UV and IR origin, these constants can be written as 


5Z g = — 4— ( -J- — -r 


2n 

SZ « = - S L 


Nf 
3 


AH^-A^j-^Ayv^), 


SZt = - 


ot s 

3tt L 


A RV (/i) - A IR (p) 
Al IV (/i) + 2A I 1 R ( / i)+4 


(2.28) 
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The contributions of the counterterms and the external self-energies to the one-loop 
matrix elements can be derived easily. For the qq channel they explicitly read 


M qq = 

J l ren 


— + — + SZ q + 5Z t ) M 
, a s m t ) 


<?<? 

LO 




2 m 


(s 35 - m 2 t ) 2 
2m t 

(s 45 -m 2 t ) 2 

( 1 


1 


+ 


, •‘>’35 - rrif S 45 — m 2 


M \f 


(2.29) 


The first term accounts for the external self-energy corrections as well as the coun¬ 
terterms of the strong and Yukawa coupling renormalization, while the remaining terms 
proportional to 8m t come from the renormalization of the top-quark mass. For the gg 
channel the former corrections have the same structure (proportional to the LO ampli¬ 
tude), but the mass renormalization involves cumbersome expressions. 

2.2.3 Pentagon diagrams 

(i) General strategy 

Four different types of pentagon diagrams occur in the calculation; a representative 
of each type is shown in the last row of Fig. [| Since the incoming quark q is massless, 
three different types of 5-point functions are generated which involve one, two, or three 
massless propagators in the loop. Each of the three types is IR divergent, so that we need 
explicit results for scalar and tensor 5-point functions in D 4 dimensions. As known for 
a long time [ |33| , in D — 4 dimensions all 5-point functions can be expressed in terms of 
4-point functions, simplifying the calculation considerably. In the following we describe 
an elegant way how this reduction can be generalized for singular 5-point functions to 
4 up to terms of order 0(D — 4) in dimensional regularization. 

I 11 order to make use of the reduction of 5-point functions to related 4-point functions, 
which works in four space-time dimensions, we first translate the dimensionally regular¬ 
ized integral E RA into another regularization scheme that is defined in D = 4 dimensions. 
To this end, we first endow all massless propagators in the loop with an infinitesimal 
mass A and call the new integral ff( mass,D ), which is identical to E ^ if A = 0. Next 
we determine a related well-defined integral, denoted E < fi™^ s,D \ with the same IR (soft 
and collinear) singularity structure as ^( mass >' D ). This integral is obtained by decomposing 
the integrand of the 5-point function in the collinear/soft limit in terms of 3-point inte¬ 
grands with regularization-scheme-independent kinematical prefactors [explicit examples 
are given below]. The difference of the two integrals, ^( mass >' D ) — Etfifi^ s,D \ has a uniquely- 
determined finite integrand and is therefore regularization-scheme independent, i.e. the 
limits ZA —> 4 and A —> 0 commute in this quantity. In total, we have generated in this 
way the relation 


E^ d) 


E ( fl + 0(D - 4) = E( mass ’ D= V - £ s ( ™ ss ’ d - 4) + 0(A). 


sing 


(2.30) 
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If artificial UV divergences are introduced via g SS,D ' ) these divergences can be controlled 
by distinguishing the space-time dimensions -Dir and D\jy for the IR and UV domains 
explicitly. As shown below, the singular parts are given in terms of 3-point functions, 
and ^( mass -- D = 4 ) can be expressed in terms of 4-point functions. Consequently, solving 
Eq. ( pR0|) for E^°\ 


E {D) = e { d) + 

^smg 1 


^gi(mass,D=4) ^gi(mass,D—4) 


sing 


+ . . . , 


(2.31) 


this integral is expressed (up to irrelevant terms indicated by the dots) in terms of 3-point 
and 4-point functions. 


(ii) Determination of singular subintegrals in 5-point functions 

Having described the general strategy, we now determine the singular terms in the 
relevant 5-point functions. The following derivation is valid for the integrals E^ mass,D \ 
which are defined for arbitrary D in the vicinity of D = 4 and which include infinitesimal 
mass regulators A for internal massless particles. The Enal formulae given for E ^ nz below 
are thus valid for both E^} and E^fg SS ' D ~ l \ We do not write the regulator A explicitly, 
since the relations for P S j ng also hold for other choices of mass regularizations^] 

We first consider the case of one vanishing internal mass, as shown in the second 
diagram of the last row in Fig. [|. for example. There is only a purely soft (logarithmic) 
divergence which is located at q —> 0 in the integrand if the vanishing mass is mo = 0. In 
this soft limit the numerator of the integral and the denominators d 2 and do are regular. 
As a result, the singularity structure of the integral remains unchanged if we set q = 0 
in the numerator and the denominators d 2 and do, since the divergence is logarithmic. 
This, in particular, shows that all such tensor integrals are finite, and only the scalar 
integral Eq develops a soft singularity, which is already contained in a scalar 3-point 
function. There are two different types of 5-point integrals with one internal massless 
line, corresponding to the second diagram of the last row in Fig. |2] and the additional 
diagram generated by shifting the Higgs line along the top-quark line. Explicitly the 
singular parts of the corresponding integrals are given by 


E 0 (~P3,Pi — P 3 ,P 4 - P 2 ,P 4 , 0 ,m t ,m t ,m t ,m t ) sing 

= 77 - 7777 - ^C o (-p 3 ,p^0,m t ,m t ), 

(ti 3 - mi){t 24 - mf) 

E 0 (~P3,Pi -P3,P4 +P5,P4,0,m t ,m t ,m t ,m t )sing 

= 77-777- 2\C 0 (-p 3 ,P4,0,mt,mt)- 

if 13 ~ mi)(s 45 - mi) 


(2.32) 


Next we consider the case with two internal massless lines, as generated, for example, 
in the third diagram of the last row in Fig. |2|. In this case, the denominators di are given 
by 


do = q 2 , d 1 = (q + p 1 ) 2 , d 2 = {q + pi - Pa] 2 ~ m 2 t: 

d 3 = {q + P 3 -P 2) 2 -m], d 4 = (q + p 3 ) 2 - m 2 t , (2.33) 

1 Alternatively, we have introduced an infinitesimal quark mass m q obeying m t rn q A > 0, and 

obtained the same final result for E^ D \ 
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where the imaginary iQ terms are suppressed in the notation. The integral becomes 
singular in the collinear limit q —> xp\ with arbitrary x, and in the soft limits q —> 0 and 
q —> —pi, which are included in the collinear case for the limits x —> 0 and x —> — 1, 
respectively. The non-singular denominators ^ 2 , 3,4 behave in the collinear limit as 

d 2 ~ (xpi +pi -p 4 ) 2 - m? = (1 + x)(it u - m?), 

d 3 ~ (xpi + p 3 - p- 2 ) 2 -m] = - x(ti 3 - m 2 t ) - xs + t 2 3 - m*, 

d 4 ~ (xpi + P 3) 2 - m 2 = - x(ti 3 - m 2 ), (2.34) 


and thus their inverse products as 

f _1_1_ (g 4 5 — t% 3) 2 _1_ 

d 2 d 3 d 4 (s 45 - mf)(fi 3 - ml) d 2 (s 45 - m 2 t ){t 13 - ~ m 2 t ){t 2 3 - m?) d 3 

1 1 
(ti 4 -m?)(t 23 -mf) d 4 

This leads to the following decomposition of the singular part of the 5-point integral, 


E...(pi,Pi ~P4,P3 ~ P 2 :P3i 0, 0, m t , m t , m ( ) sing 

= 7 -2577- ^C.SPhPl - P 4 , 0, 0, 77lt) 

(s 45 -mi)(t l3 - ml) 

_ (s 4 5 ~ t 2 z) 2 _ 

(s 45 - ml)(ti 3 - m 2 ){tu - ml)(t 2 3 ~ ml) 

+ (t u -™;43-™?) c - (Pl ' P3 ' 0 ' 0 ' m,) ' 


C...(pi,Pa~P2,0,0,mt) 


(2.36) 


This relation is valid for all tensors; in contrast to the case with only one vanishing 
mass, in this case also tensor integrals become singular. Note that E sing does not only 
include the singularities but also regular terms in a well-defined way, which is important 
when switching from one regularization scheme to the other. 

Finally, we consider the case with three internal masses zero, as encountered, e.g., in 
the first and last diagram of the last row in Fig. |j. Here three soft and two collinear limits 
give rise to divergences. Identifying m 0 = m 4 = m 4 = 0, the collinear limits are realized 
by setting q —> xp\ or q —> yp 2l and the three soft limits are reached for particular values 
of x and/or y (i.e. x = y = 0, x = —1, or y = 1). The decomposition for the two collinear 
limits can be worked out in analogy to the previous case. When combining the two limits, 
however, double-counting of the soft singularity at x = y = 0 has to be avoided. The 
decomposition of the singular structure of the 5-point integral finally reads 


E,„(p 1 ,p 1 -p 3 ,P4 ~p 2 , —p 2 , 0, 0, m t , m t , 0) sing 

= 77 m 2 )(t ztit — ^ 2 ,0,0,0) 

[t 13 — m t ){t 24 — m t ) 

+ T7--- —C...(pi,Pi — P3, 0 , 0 , m t ) 


s ( s 35 - ml) 


(t'24 ~ ^ 35 )^ 


s(s 35 - ml)(ti 3 - m 2 t ){t 24 - ml) 


C...(pi,p 4 


p 2 , 0, 0, m t ) 
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1 


+ 


s(s 45 - mf) 


('- Ih- ~P 2 , 0, m t , 0) 


(tl 3 — 545)^ 


s(s 45 - rn 2 )(t 13 - m 2 )(t 2 4 - m 2 ) 
which is again valid for all tensors. 


C'...(Pi-P3,-P2,0,mi,0), (2.37) 


(to/ The scalar 5-point function 

In 77 = 4 space-time dimensions the scalar 5-point function E 0 can be expressed [ [33|] 
in terms of the five related 4-point functions D 0 (i), obtained from E 0 by omitting the 
denominator <7*, as defined in Eq. (|2.24|). Using the shorthand notations 


Wo = Y i0 = Y 0i = m 2 0 + m 2 - q 2 , Y %J = m 2 + m 2 - {q, : - qj) 2 , i, j = 1, 2, 3,4, 

(2.38) 

the relation reads 



7 -! (mass, D=4) 

“ A) 

4"4«) D<">(1) B<->(2) 

44 3 ) 

444 ) 


1 

boo 

Wi 

W 2 

w 3 

w 4 

0 = 

1 

b '10 

Wi 

W 2 

w 3 

w 4 


1 

Wo 

Wi 

W 2 

w 3 

w 4 


1 

Wo 

Wi 

W 2 

w 3 

w 4 


1 

Wo 

Wi 

W 2 

w 3 

w 4 

where the scalar integrals are defined for D 

= 4. Solving for Eq 

nass,D=4) 


7 -j(mass ,D= 

=4) _ 

1 

£ det(W)^ mass ’ D= 

4 = 0 

=4) ru 


-^o 


det(W) 



(2.39) 


(2.40) 


where Y = (Y/,), and Y t is obtained from Y by replacing all entries in the it\i column 
by 1. 

Although Eq. (|2.39|) was derived for D = 4, it is valid (in all cases considered here) 
even for arbitrary D up to terms of order 0(77 — 4): 


d D| = - £ det(y f ) + 0(0 - 4). (2.41) 

This generalized relation simply follows from the result derived in \( ii )\ that the same linear 
relation between E s \ ng and 3-point subintegrals holds in all mass/dimensional regulariza¬ 
tion schemes. 


(iv) Tensor 5-point functions 

We evaluate all tensor integrals by calculating the coefficients to the Lorentz covariants 
that span the tensor. For the pentagons, tensor integrals E^ v appear up to rank 4. In 
practice, however, these integrals are needed only up to rank 3, since one momentum 
factor q in the numerator of a top-quark propagator can always be cancelled against 
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a propagator denominator. The covariant decomposition of the relevant 5-point tensor 
functions reads 


y ] Qi,p.Ei, 


i= 1 
4 

Ep,v y ) Qi,pQj,uEij T Qfj,vE 00; 

i,j =1 

4 4 

E'fj.vp y ) Qi,pQj,vQk,pEijk + E ^Sti,p9vp T Qi,u9pp T QiypQpLv) ^OOi- 

i,j,k= 1 i=l 


(2.42) 


Since the tensors E^^ are symmetric, the tensor coefficients E tJ and E ljk are symmetric 
in the indices. The coefficients E i: E i3 , etc., can be recursively reduced to the scalar inte¬ 
grals Eq, Do, etc., by an algorithm known as Passarino-Veltman reduction, as described 
in Refs. [ 0, or by special techniques [ |34|, [35j developed for higher n-point tensor 

integrals. The resulting expressions are cumbersome and do not reveal any new insight; 
therefore we do not list the expressions explicitly. 

Note that the metric tensor appearing in Eq. (|2.42|) is D-dimensional. In contrast to 
the D-dimensional metric tensor g a g , its four-dimensional part, g^g~ 4 \ is spanned by the 
linearly independent base vectors qf 


(-0=4) 
9 a? 


E { z % 


z ij = ( QiQj ), 


(2.43) 


i,j =1 


where Z is the Gram matrix of the momenta g*. The difference A g a p = g a g — ^ 

thus, of 0(D — 4). It is interesting to notice that the terms proportional to A g a g do not 
contribute in our case. Since A g a g is an object of 0(D — 4) it could only contribute if an 
accompanying tensor coefficient, such as E 00 , is divergent. The explicit tensor reduction 
to scalar integrals, however, reveals that all such coefficients are finite. This fact can 
also be derived from the following argument. In tensor integrals the only covariants that 
come with divergent coefficients are built up by the singular regions in momentum space. 
Therefore, a collinear divergence occurring for q —> xp can only show up in tensors built up 
by the momentum p alone. For purely soft (logarithmic) divergences occurring at q —> 0 
power-counting shows that only the scalar integral is divergent; shifting the integration 
momentum by p, i.e. q —> q + p , only leads to divergent coefficients related to tensors of 
p alone. UV divergences could occur in covariants containing A g a p] however, all relevant 
5-point integrals are UV-hnite, which is easily proven by power-counting. 


(v) Alternatives in practice 

The results on 5-point functions described above admit two conceptually different ways 
to evaluate the pentagon diagrams. We followed both ways, as described below, and we 
found perfect numerical agreement. 

In the first approach the scalar integral Eq is evaluated in D dimensions using Eq. (p.41|) 
and also the reduction of the tensor coefficients defined in Eq. (|2.42| ) is performed consis¬ 
tently in D ^ 4. In this approach all divergent scalar 3- and 4-point integrals Co and D 0 
are needed in dimensional regularization. We have listed them in Appendix [A]. 
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Alternatively one can make consistent use of Eq. ([2.31 ) and its analogous version for 
4-point functions, in order to replace all D-dimensional 4- and 5-point integrals by four¬ 
dimensional integrals and D-dimensional 3-point functions that contain the IR (soft and 
collincar) singularities. The reduction of the 4- and 5-point tensor coefficients then works 
in four IR dimension^. An advantage of this procedure is the easy analytical control of 
all IR singularities which are contained in 3-point functions only. In this approach all 
divergent D 0 functions are needed in the mass-regularization scheme, and all divergent 
Co functions in both the dimensional and the mass-regularization schemes. The mass- 
regularized functions can be derived from the results of Appendix |AJ Another benefit of 
evaluating pentagon diagrams in four dimensions is the direct application of the tensor 
decomposition of Ref. [|35|. This procedure consistently avoids leading Gram determi¬ 
nants (i.e. those of four momenta) that are potential sources of numerical instabilities, as 
described in the next section. 


2.2.4 Numerical stabilization 

As explained above, all tensor integrals occurring in the calculation can be algebraically 
reduced to scalar integrals following the well-known Passarino-Veltman algorithm [ [28|] . 
This procedure is based on a decomposition of the tensor integral into a complete set of 
covariants that are spanned by the metric tensor and all independent momenta involved 
in the integral. The coefficients of the covariants are obtained by inverting a set of linear 
equations, which is non-singular as long as the momenta qi spanning the covariants are 
linearly independent, i.e. as long as their Gram determinant G = det(g^) is non-zero. 
Each tensor rank adds a factor G _1 to the expression of the tensor integral. Near the 
boundary of phase space, some determinants G become arbitrarily small, since some mo¬ 
menta become linearly dependent at the boundary. The huge factors G~ n are balanced 
by cancellations between different tensor coefficients that are calculated in terms of scalar 
integrals; the whole virtual correction remains well behaved at the phase-space boundary. 
These cancellations lead to severe numerical instabilities in the evaluation of the correc¬ 
tion. The instabilities are most pronounced in the pentagon diagrams, which do not only 
possess the highest tensor ranks, but also involve additional sources of cancellations. For 
instance, the decomposition ( |2.43| ) of the four-dimensional metric tensor also leads to a 
Gram determinant in the denominator, and the reduction (|2.39| ) of the scalar five-point 
function Eq in terms of five four-point functions involves Caley determinants det(Y)..) that 
can become very small. 

If the virtual correction is not evaluated with particular care at the phase-space bound¬ 
ary, a phase-space point will soon be selected in the Monte Carlo integration where the 
result for the correction is totally wrong. The wrong result will typically be a large 
number owing to incomplete cancellations in the evaluation. Such unphysically large con¬ 
tributions occur more and more frequently in the adaptive integration procedure and they 
thus destroy the result of the phase-space integration completely . 

When using the Passarino-Veltman procedure, we avoid the numerical instabilities by 
extrapolating the virtual correction from the numerically safe inner region of phase space 

2 UV divergences, which artificially occur in subexpressions during the tensor reduction, still have to 
be treated in D dimensions. 
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to the numerically delicate boundary. To this end, we first divide the full phase space 
into two regions, one where the correction is evaluated directly and another where the 
correction is obtained from an extrapolation out of the first region. 

The extrapolation procedure is applied to the full one-loop correction, since any sub¬ 
set of diagrams (such as the pentagons alone) in general shows stronger variations over 
the phase space than the full correction, owing to cancellations between different dia¬ 
gram types. Therefore, we have to inspect the whole boundary of phase space, which is 
characterized by the zero of the Gram determinant G of the four momenta pi ,... ,p 4 , 
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(2.44) 


where k is defined in Eq. (|2.21|) and p§ is fixed by momentum conservation, p\ + P 2 = 
p 3 + Pi + P 5 - The three-momenta p, and the angles 9 1 , y are defined in the partonic 
centre-of-mass frame. The angle 9± is the angle between P 3 and the beam axis, and y 
denotes the (azimuthal) angle between the p 4 ps plane and the plane spanned by P 3 and 
the beam axis. We identify the “safe” region by demanding G/s 4 as well as siny to be 
larger than certain numerical cut values. The extrapolation into the “problematic” region 
is done as follows: 


1. Regarding the integrand as a function of the variable y, thereby keeping p 3 and 
| p x | fixed, we expand it in terms of a Fourier series. As a set of basis functions we 
took trigonometric functions or Legendre polynomials. The first few coefficients of 
the series are determined by numerical integration, for which we took Simpson or 
Gaussian quadrature. If the extrapolation result is close to the directly calculated 
result for the correction, the direct result is taken over. 

2. If the agreement between direct and extrapolated result is too bad, the first step 
is repeated again with twice as many coefficients in the Fourier series and twice as 
many points in the related integration for the coefficients. The first step is repeated 
at most twice. If the directly calculated result is still not in good agreement with 
the extrapolation, but instead the various extrapolations look convergent, the most 
precise extrapolated result is taken over. 

3. If the procedure still has not stopped yet, step one is modified by changing the 
extrapolation direction by varying also p 3 and |p 4 | during the change in y. Step 
two is carried out as above. 

4. The modified extrapolation of step three is repeated if still no success has been 
reached. Finally, if none of the extrapolations looks convincing, the integration 
point is discarded, but this does not happen very often. 

The extrapolation procedure solves the problem of numerical instabilities at the expense 
of computing time. Of course, the described procedure has to be optimized by selecting 
appropriate cut values on G/s 4 and siny and by choosing appropriate techniques for the 
calculation of the Fourier coefficients. Nevertheless the final phase-space integration has 
to be stable against any change in the extrapolation procedure. We have checked this 
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stability by changing the cut parameters and by choosing different systems of orthogonal 
functions as well as different integration techniques, as mentioned above. 

In addition to adopting the usual Passarino-Veltman reduction of 5-point tensor co¬ 


efficients, we have applied the method of Ref. |35| where the strategy [ |33| for scalar 
5-point functions is generalized and 5-point tensor coefficients are reduced to 4-point inte¬ 
grals directly. In this approach inverse Gram determinants of four momenta are avoided 
completely. Applying this alternative renders the virtual correction near the phase-space 
boundary much more stable than in the usual Passarino-Veltman approach. The results 
obtained by the two methods mutually agree with each other. 


2.2.5 One-loop soft and collinear divergences 

There are two sources of IR (soft and collinear) singularities in the one-loop diagrams: 
the self-energy corrections of the external fields and various kinds of 3-point subintegrals 
whenever a (nearly on-shell) gluon is exchanged between external lines. The former 


singularities can be read off from the results of Section |2.2.2| , and the strategy to determine 
the latter singularities has already been explained in Section |2.2.3| \(v\ . Here we list the 
explicit results: 


MV- loop 

M{ 9 _ 


loop 


IR 


IR 


Qf s 

247T 


(/?? - T/fl - 2/|| - 2/« - 7/ff + /??) Mfo, 


(2.45) 


cr s 

247T 


(9/ff + 9/ff + 9/ff — /|f) ( -Cf + Cf + Cf ^ ^AT direct + -M fusion 


+ (9/ff + 9/ff + 9/ff - /If) (J<T - Cf + Cf) (M crossed - l -M iusion 


+ 6(/ff - fll - /f| + /|f) Cf (M direct + -M fusion 
+ 6(/if - /if - /If + /If) Cf ^M crossed - ^Affusion 


(2.46) 


The functions AT... are given in Eq. (|2.19|) , and the divergent factors ,/f and ff? read 

3' 


/ff = 2Af(p) - 2Af(p) 
/f| = 2A®(/x) - 2A I 1 R (/i) 


In-77 — iO — 

V ) 2 


, / s . \ 11 AT 

l m\ ) 6 9 


m = a® (p) - a ir (p) 


2 In 1 


tjj 


Jif = Af (/r) - Af(p) 

fi! = fil = 2AF(/0 


A; 


2 In I 1 - 

V 


_ 17 A^ 

~ y + t 


s 34 - 2 m 2 
5 34/3 


S34 


4 ffif ^ 34 ] ) + 1 

\As34 + 1/ 


* = 1,2, / = 3,4, 


A*34 — \/ 1 


4m| 

S34 


+ iO. (2.47) 


These results are in agreement with Ref. [0, where the general singularity structure of 
QCD amplitudes has been presented. The agreement becomes apparent after rewriting 
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Figure 3: Generic set of Feynman diagrams contributing to real gluon emission and gluon- 
(anti)quark scattering. 


Eqs. (|2.45|) and (|2.46|) in the following generic form which is valid for both the qq and gg 
channels, 


M !_ 


loop 


ry 

m = M l0 <8 £ (T.. T,) /„, 


(2.48) 


i,j =1 
i<j 


36|1 , and the symbol ® denotes 


Here T.,. are the colour operators defined in Refs. [ ^2 
the colour correlations. Since we will make repeated use of this notation in the remainder 
of our report, a short description can be found in Appendix |B]. 


2.3 Real corrections 

The real corrections at relative order 0(a s ) include processes with real-gluon radiation 
as well as ttH production in gluon-(anti)quark scattering. After a brief description of the 
parton reactions and their evaluation, we will discuss in some detail the isolation of 1R 
divergences for soft and collincar parton configurations as well as the treatment of collincar 
initial-state singularities by means of mass factorization. 


2.3.1 Parton processes and their evaluation 

The evaluation of the 0 (a^ g^ tH ) cross section requires the calculation of the gluon 
bremsstrahlung processes 

?(Pi) + q(P2) tiPs) + t{p±) + H(p 5 ) + g(p 6 ), 

g(pi) + g{P2) -> t(p 3 ) + t{p 4 ) + H(p 5 ) + g(p 6 ), (2.49) 

and the gluon-(anti)quark scattering reactions 

q(Pi) + g{P2) -> t{p 3 ) + t(p 4 ) + H(p 5 ) + q(p 6 ), 

g(pi) + ?(P 2 ) -> t{p 3 ) + i(p 4 ) + H(p 5 ) + q{p s ). (2.50) 


A representative set of Feynman diagrams, contributing to the real corrections, is given 

in Fig. [i. 

Hclicity amplitudes for the processes (|2.49| , p.50| ) have been generated and evaluated 
using the program packages MadGraph [ 37] and HELAS [ [38|. The result has been 


verified by an independent calculation based on standard trace techniques. 
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2.3.2 Soft and collinear divergences 


The 2^4 cross sections contain singularities in the limits where the splitting g —» 
gg, 9 QQ ; 9 —> 99i or q —> qg becomes soft or collinear. We have used the dipole 


subtraction method [|23|], as formulated in Ref. [ [24|| for massive QCD partons, to extract 
the singular part of the real cross section in D dimensions. The dipole subtraction method 
is based on the fact that soft and collinear limits of a real-emission matrix element can 
be expressed by convoluting process-independent splitting functions with the LO matrix 
element. By subtracting these splitting-function expressions as counterterms from the 
squared real matrix element, the 1R (soft and collinear) singularities cancel and the phase- 
space integration can be performed numerically in four dimensions. The counterterms 
are constructed such that they can be integrated analytically in D dimensions over the 
phase space of the extra emitted parton, leading to poles in e = (4 — D)j 2. When the 
integrated subtraction counterterms and the virtual cross sections are combined, the soft 
singularities cancel. The remaining initial-state collinear divergences have to be absorbed 
into a redefinition of the parton distribution functions at NLO. Below we will describe 
the subtraction formalism in more detail and we will give explicit results for the soft and 
collinear contributions. 

The NLO contribution of the cross section can symbolically be written as [ 
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0nlo(Pi»P2,/4’) = <7nlo } (Pi>P2) + (Tnlo } (Pi ; P2) 


dx 


d^rf(-Pi, ■ • • ,Pe) 


~ab{ 3}/ 2\ i ~ab{ 3}/ 2\ 

%o (h x PliP2, Pf) + <Wo ( x ; Pl» 9f) 


Je=0 


E( dcj : 

L dipoles 


ab 

LO 


dhdipole) (pi, ... ,P5,Pe) 


J e=0 


+ 


' 3 L 


d (Ty (pi, • • • ,Ps) + do£ 0 (pi,... ,p 5 ) ® T 


rab 


E 

'=9.5 

E 


da: 


“'=9.9,5 0 


/ 3 L 


da: 


=9.9.9 


/ 3 L 


d^Lo(zpi,P2,P3,P4,Ps) 
d ^Lo(Pl^P2,P3,P4,P5) 1 


Je=0 


ttH-\-b 


x , ii 2 f ) 


Je=0 


T 


b,b' 

ttH-\-a V 


( x ,P 2 f) n • (2.51) 


Je=0 


This symbolic expression is valid under the assumption that no massless particles are 
identified in the final state, i.e. the additional gluon or light (anti)quark is integrated out. 
The indices a and b denote the flavours of the incoming partons, i.e. ab = qq, gg, qg , gq in 
the processes (|2.49| , pT5C|) , and hf is the factorization scale at which the redefinition of the 
parton distribution functions is performed at NLO. The notation for the integrals indicates 
that the real-emission contribution (labelled by R) involves four final-state particles [one 
gluon or light (anti)quark more than in LO], while the virtual contribution (labelled by 
V) has LO three-particle kinematics. 

The dipole subtraction counterterm is given in terms of universal splitting functions, 
dVdip 0 ie (see Ref. [ ^4|] for explicit expressions), and an appropriate colour, spin, and flavour 
projection of the LO cross section, dcrLO- The symbol eg) denotes such possible colour, spin, 
and flavour correlations. In particular, in the sum over “dipoles” the flavour correlations 
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can change the incoming partons a or b in dcr£o- The final state of da^o is fixed to ttH , 
since no light (anti)quark can be radiated from the top quarks. As far as the four-particle 
phase space is concerned, each individual dipole counterterm is evaluated for a specific 
set of auxiliary momenta pi,... ,ps. This involves combining the emitting particle (with 
momentum pj) and the emitted light particle (with momentum p g) into a single parent 
particle called “emitter” (with momentum pj). In order to balance the total momentum, 
one of the other external particles (with momentum p/,) is replaced by a “spectator” 
particle (with the same quantum numbers, but modified momentum pf.). Apart from the 
situation where both emitter and spectator are incoming partons, the momenta of the 
remaining external particles are not altered]]. This phase-space mapping is constructed 
in such a way that 

- the soft and (quasi-)collinear momentum configurations are reproduced point-wise; 

- both the emitter and spectator are on their respective mass shells; 

- momentum conservation is preserved; 

- the phase space can be factorized into a hard-scattering phase space and a single- 
particle phase-space factor corresponding to p e . 


Many of these phase-space mappings either have an initial-state emitter or an initial- 
state spectator. In those cases the corresponding initial-state momentum is rescaled by 
a longitudinal momentum fraction x, resulting in a boosted reference frame for the final- 
state momenta and a phase-space factorization that involves a convolution over the boost 

1231, 


m 


parameter x. For explicit phase-space mappings we refer to the literature 

The readded subtraction counterterm, integrated over the phase space of the emitted 
extra parton, can be split into two parts. The first part, involving the universal insertion 
operator X ab , has LO kinematics and contains all the poles in e that are necessary to cancel 
the soft singularities in day . These so-called endpoint contributions are obtained by 
integrating the dipoles over the complete phase space of the emitted extra parton and by 
using the (•••)+ prescription to extract the singular terms related to the soft endpoint (x —> 
1) of the x-integration of the dipoles. The second part is a finite remainder that is left after 
factorization of initial-state collinear singularities into the parton distribution functions at 
NLO. It involves LO cross sections with boosted three-particle kinematics, x-dependent 
structure functions T(x,p,p), and an additional one-dimensional integration with respect 
to the longitudinal momentum fraction x of an incoming parton. The structure functions 
T(x,/4) depend on the factorization scheme; only the terms proportional to In (p/p) are 
factorization-scheme independent. The summations over the flavours a' and V represent 
the afore-mentioned flavour correlations. Non-zero contributions are obtained if both 


d<T£o and X^ a H+b are non-zero, or if both da£o and V b t f H 


ttH-\-a 


are non-zero. 


In the remainder of this subsection we give the explicit expressions for the readded 
subtraction counterterms, which are based on the results of Ref. [24], The explicit con¬ 


struction of the dipole terms, which is not described here, can also be found there. We 


3 In the case of both an initial-state emitter and an initial-state spectator, the momentum of the 
spectator is left unaltered. Consequently, the momenta of all final-state particles undergo a Lorentz 
transformation. 
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start with the insertion operator X ab , which depends on the colour charges, momenta, and 
masses of the final-state particles in the LO matrix element: 


X” = X" = 0, I”* = - £ E (Ti ■ T,) J‘ b , ab = qq, gg, 


(2.52) 


i,j =1 
i<j 
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(2.53) 


The variable /3 S34 is defined in Eq. ( [2.47]) and p is given by p = (1 — As 34 )/( 1 + /3 S34 ). The 
dilogarithm Li 2 (x) is defined in the usual way as 


Lio(x) = 


rl ln(l — art) 


jo t 

The Casimir operators T 2 and anomalous dimensions y a read explicitly 

4 X, "In 3 


(2.54) 
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f 
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(2.55) 


Inserting all this in Eq. (|2.53|) , we observe indeed that the singular Al, R (/i) and A* R (/i) 
terms in the functions cancel against the singular terms in the coefficients /“ b given 
in Eq. (|2.47p . 

Finally, we give a general expression for the structure function r// +fe (x, /i/) in the 
MS factorization scheme: 


r tw+b = In ( 4 ) P aa \*) + 2P r a e “'(x) ln(l - x) + P' aa ,{x) + 4<5““'T 


aa' rr\ 2 
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Pd 


4(1 - x + /x-i) 2 _ 

(2.56) 


The various Altarelli-Parisi splitting functions P aa ' (x), P r “g(x), and P' aa i{x) can be found 
in Appendix [C], and the kinematical variable p 2 ]X is given by 


Pd = 


mr 


2piPj 


(2.57) 


where pj is the auxiliary momentum of particle j in the boosted reference frame, while pi 
is the original initial-state momentum of the incoming parton a. In our convention the 
(...) + prescription is defined by 



d*h(x) /(.r. //;) ! o($(x)) 


dx / d<f>(x) / dp, 6{fi-fj$) f(x,p) g($(x)) 


= J o dx J^dpf{x,p) !^J d$(x)gmx))6(p-pff) 

-/d4(l) S (4(l))^-g 1 1 , )| 

- <1.7. | (d'Kx j /(l./i'fi i - /d'lM'l ) Jix.lpn flCI’M j j) , (2.58) 


where the label (x) of p!'^ indicates that pj belongs to the phase space <f>(x) of the x- 

boosted reference frame and therefore p^ = m t /\/‘2pipj. In fact, this is the same as saying 
that pji should be kept fixed during the x-integration, e.g. by using it as a phase-space 
integration variable in d<3>(x). The expression for the structure function T^ +a (x, p 2 F ) can 
be obtained from Eqs. ([2.56D and ( |2.57| ) by interchanging (a, a 1 , pji,pi) and (b,b', p,j 2 ,P 2 )- 
The formulas presented above are spelled out for total cross sections, but they can 
also be used to evaluate any kind of IR-safe differential distributions by filling histograms 
during the Monte Carlo integration. In this context, the only subtlety in the dipole sub¬ 
traction approach is due to the fact that contributions of different dipole terms in general 
contribute to different histogram bins. The IR safety of an observable guarantees that 
singular terms of the original integrand and their corresponding counterterms contribute 
to the same bin in the soft and collinear limits.^ 

4 If a direct calculation of a differential cross section is aimed at, i.e. without an integration over the 
full phase space, further effort is needed, since the kinematics in some dipole terms would have to be 
changed for this purpose. 
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2.3.3 Slicing method for the qq channel 

In the phase-space slicing approach (see Ref. [ [39j for a review) the soft and collincar 
regions are excluded from phase space by appropriate phase-space cuts, generically de¬ 
noted by A cut . The singular integration over these regions is then carried out analytically 
using factorization properties. The full real correction is obtained by the sum of these 
two contributions in the limit A cut —> 0. In practice, a plateau is required for small but 
finite values of A cut . 

In the following paragraphs we describe the so-called two-cutoff method for the qq 
channel, where the soft region for gluon bremsstrahlung is defined by a cut A E on the 
gluon energy and the collincar region for gluon radiation off the initial state by a cut A 6 
on the emission angle. 

(i) Soft-gluon region 

Denoting the momentum of the radiated gluon by k, the soft region is defined by 


0 < k 0 < A E -C m t , 


(2.59) 


where k$ is the gluon energy in the partonic centre-of-mass frame. In the soft region, gluon 
radiation is described by an eikonal current (see e.g. Ref. [ Q), and the contribution to 
the bremsstrahlung cross section is given by 

a 4 

d^soft = - cKo ( T * • T i) 9ij(Pi,Pj), (2.60) 
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where gij denote the soft integrals 
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The integrals are obtained as follows: gi 2 is calculated easily; g 34 can be taken over from 
Ref. [ [3l]], where this integral is given for an infinitesimal gluon mass m g that translates 
into 1/e via the substitution ln(m 2 ) —*• (47r/r 2 ) e T(l + e)/e; the remaining g tJ can be 

derived using the auxiliary integrals (C.20) and (C.25) of Ref. [ |4(J. The results are 
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The 1/e 2 poles, i.e. the A!, R (/i) terms, of the functions g l3 are the same as in the coefficients 
fij given in Eq. (|2.47|) . This implies that these terms cancel in the sum of virtual and 
soft corrections, as it should be. 

(ii) Collinear gluon emission from the initial state 

The region of collinear gluon emission from the initial state is defined by 

k 0 > A E, 0 < 0(p t , k) < Ad < 1, i = 1,2, (2.65) 

in the partonic centre-of-mass frame. According to the QCD factorization formula (see 
e.g. Ref. [ 23|), the contribution of collinear gluon radiation off parton i with i — 1,2 
reads 
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( 2 . 66 ) 


where the Casimir operator C-p = 4/3 has been defined in Eq. (|2.55|) and the Altarclli- 
Parisi splitting function P qq (x ) can be found in Appendix 0 These results can, for 
instance, be easily derived by evaluating the x and v integrals in Section 5.5 of Ref. [ 
with the A E and A 9 cuts. 

From the results given above it is easy to verify that all divergences Af R (/r) cancel i 


m 


the sum of the virtual corrections (see Section |2.2.5|) , the soft corrections dof^, and the 
part of the collinear corrections d<rff u with the LO kinematics, i.e. the part proportional 
to dcrfoGoi). The remaining P qq (x) Af R (/r) term is absorbed into the redefined parton 
densities. 


(in) Independent check within the mass-regularization scheme 

The sum of dof,f ft + dcrff u has been checked by deriving all formulae for infinitesimal 
quark masses. Apart from g 34 this calculation is completely independent to the procedure 
in paragraphs (i) and (ii). These are the steps to be followed: 


1. Derive derf^ and daf^ for infinitesimal quark masses for q and q ; see Ref. [ [H|] for 
41j, in particular Eq. (26), for dcrff u . 


dofo ft , and Ref. 


2. The translation of the continuum part of duffn [the part with P qq (x) and the (...)+ 
prescription] to the dimensional regularization scheme can be performed by com¬ 
paring the renormalization of the parton distribution functions within the two reg¬ 
ularization schemes. Specifically, Eq. (6.6) of Ref. [ has to be compared with 
Eq. (23) of Ref. [ |42j (see also Ref. [ |43|). 


3. The translation of the endpoint part of dof^ + dcr^ u [this is the part proportional 
to daLo(pi)] is described in Eqs. (23) and (41) of Ref. 


2.4 Leading threshold corrections 

In view of the complexity of the NLO QCD corrections, as detailed in the previous 
subsections, any additional cross-check is welcome. The partonic energy regime that 
is most suited for performing such cross-checks is the threshold region, \/§ > 2 m t + 
Mh , where the QCD corrections are dominated by soft-gluon effects. The leading soft- 
gluon effects fall into two categories: (virtual) Coulombic gluon exchange between the 
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Q,9 


q,9' 

Figure 4: Generic diagram leading to the Coulomb singularity 



massive final-state particles and IR-enhanced edge-of-phase-space effects from (real) gluon 
radiation. Both effects can be calculated analytically. 

Apart from providing explicit cross-checks, the leading threshold corrections can also 
serve as a simple, intuitive tool for understanding the differences between the hadronic 
NLO corrections at the Tevatron and the LHC. This should not be overrated, though. In 
ttH production, with its massive three-particle phase space, the influence of the threshold 
region is substantially less than for processes with massive two-particle final states, like 
tt production (see Section |3T|) . 


2.4.1 The Sommerfeld rescattering correction 

For energies near the threshold the virtual QCD corrections are enhanced by Coulorn- 
bic gluon exchange between the top and antitop quark in the final state. The corre¬ 
sponding generic diagram is shown in Fig. f|. When the gluon momentum tends to zero a 
so-called Coulomb singularity is generated. In this limit the contribution of the singularity 
to the matrix element, denoted by Alcouh is given by 

jrry 

M Coul = - Mlo ® ttt-~(T 3 • T 4 ), (2.67) 

4ps 34 

and it increases inversely proportional to the t/t velocity (3 S34 in the tt centre-of-mass 
frame, as defined in Eq. ( j2.47]) . This is the well-known Sommerfeld rescattering correction 
. Upon averaging 1 / (3 S34 over the phase space near threshold, we obtain the following 
singular contribution to the partonic cross section 
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( 2 . 68 ) 
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The coefficient (3 is proportional to the maximal relative velocity of the tt system. The 
rescattering coefficient Cc ou i depends on the production mechanism of the tt pair, 
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(2.70) 


In e + e~ annihilation the tt pair is generated in a colour-singlet state in which the quark 
and antiquark attract each other, and with = +4/3 the correction is positive. This 

leads to a strong enhancement of the e + e _ —> ttH annihilation cross section near threshold 
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By contrast, in qq annihilation the tt pair is generated in a colour-octet state in 
which the force is repulsive, and with Cq ou1 = — 1/6 the correction is negative and 
relatively small. This configuration is realized in the dominant channel at the Tevatron, 
qq —> ttH. Finally, in gg fusion the various colour channels contribute differently, and the 
relative weighting depends on the ratio M^/Am^. For not too large Mh the rescattering 
coefficient Cq ou1 is moderately positive, and for Mfj/Arn^ <C 1 the rescattering coefficient 
approaches Cq ou1 (M h = 0) = 11/42 for gg —■> tt (see Refs. [ P>f). 

In the close vicinity of threshold, our numerical results for the NLO partonic cross 
sections indeed reproduce these Coulomb singularities (see Section |XT|). 


2.4.2 Threshold logarithms from soft-gluon radiation 

Owing to the lack of available phase space near threshold, enhanced real-radiation 
effects from the initial state (ISR) and the final state (FSR) are only noticeable in the 
vicinity of soft IR poles. In that case the phase-space suppression is compensated by 
the inverse dependence on the energy of the emitted massless particle, resulting in an 
integrated correction that is logarithmic in the phase-space volume - in contrast to the 
Rutherford pole in the Sommerfcld correction. The calculation of these so-called threshold 
logarithms simply amounts to a convolution of the LO cross section with the soft-pole 
parts [oc 1/(1 — a:)] of the dipole functions for real-gluon emission. 

The LO matrix elements of the reactions qq, gg —> ttH are basically constant near 
threshold. This leads to the following threshold behaviour of the LO cross section [using 
Eqs. (127201) and (gjjp ]: 
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where W^q represents the part of the LO cross section that remains constant near thresh¬ 
old. The strong suppression of the LO cross section near threshold oc /3 4 is a consequence 
of the massive three-particle phase space. 

The threshold logarithms at NLO are given by 
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The function Fi og (x) comprises all IR-sensitive parts of the structure functions 
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where ///? denotes the factorization scale. Exploiting charge conservation and keeping only 
the terms involving In (3, one obtains finally 
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and 


Qsr — 2 C F — -, Cf s 9 R — 2 Ca — 6, 
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(2.76) 


The contributions from final-state gluon radiation exhibit a behaviour in colour space 
that is parallel to the Coulomb singularities. In qq annihilation the tt pair is generated 
in a colour-octet state, resulting in an overall colour charge Cf| R = 3. In the case of gg 
fusion both colour-octet and colour-singlet states are superimposed. However, soft gluons 
cannot be emitted from colour-singlet tt final states near threshold owing to the overall 
vanishing colour charge and the long wavelength of the emitted IR gluons - as evident 
from e + e - annihilation [[^J. Therefore only the colour-octet state receives a logarithmic 
correction, causing the overall coefficient Cpf R to depend on the ratio M^/4m 2 . 

In total, the threshold logarithms give rise to a positive NLO correction near threshold, 
as long as the factorization scale is not chosen too small. Furthermore we observe that the 
corrections in the gg channel are larger than the corrections in the qq channel by roughly 
a factor of two. A threshold analysis of our numerical results for the partonic NLO cross 


sections indeed confirms the above-given threshold logarithms (see Section |3T|) . 


3 Numerical Results 


In this section we present the NLO QCD results for the total partonic cross sections 
and for the final hadronic cross sections at the Tevatron and at the LHC, including the 
differential distributions in transverse momentum and rapidity of the Higgs boson and 
of the top quarks. As discussed previously, renormalization and factorization have been 
performed in the MS scheme with the top-quark mass defined on-shell. The top quark 
is decoupled from the running of the strong coupling a s (n) and from the evolution of 
the parton densities. For the evaluation of the hadronic cross sections we have adopted 
the MRST [ [47| parton densities at LO and NLO, corresponding to the QCD parameters 
Ag 0 = 167 MeV and A^ s = 239 MeV at the one- and two-loop level of a s (p), respectively. 
The strength of the SM Yukawa coupling is fixed by g tt H = mt/v, where v = (v^G^) -1 / 2 
is the vacuum-expectation value of the Higgs field and Gf = 1.16639 x 10 -5 GeV -2 . The 
top-quark mass is set to m t = 174 GeV. 


3.1 Partonic cross sections 

To explore the NLO QCD results at the parton level for the production of ttH in 
quark and gluon collisions, we introduce the scaling functions f a b{j])i 


= //(«(,,) +4TO , ( y) 
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The centre-of-mass energy of the partonic reaction, \/I, is expressed in terms of the 
quantity 77 = s/ 4/Xq — 1 , which is better suited for analyzing the scaling functions in the 
various regions of interest. The reference scale p 0 — + A/#/2 has been set to half the 

partonic threshold energy. The index pair a, b indicates the partonic q, q , g initial state 
of the reaction ab —> ttH + A". We have identified the renormalization scale with the 
factorization scale: /Pr = Pf = / 1 . The scaling functions are divided into the LO term 
the scale-independent NLO term C 1 ) anc } the scale-dependent NLO contribution 

/ (1) MTV/dj)- 

The scaling functions are displayed in Fig. ^]for the quark-antiquark, gluon-gluon, and 
(anti)quark-gluon channels. The Higgs-boson mass is set to M H = 120 GeV. As described 
in Section |2.4j , the parton cross section near the production threshold (77 ~ (3 2 -C 1) is 
dominated by soft-gluon effects which can be calculated analytically. The inserts in Fig. 
show the ratio of the threshold approximation, cr t hres = crcoui + ^iog, Eqs. ( |2.68|) and ( p.74|) , 
and the numerically integrated, complete NLO scaling functions. It is evident that the 
complete NLO calculation approaches the threshold approximation as 77 —^ 0. However, 
because of the (3 4 suppression of the massive three-particle phase space, see Eq. ( |2.71|) , the 
impact of the Coulomb singularity and the threshold logarithms from soft-gluon radiation 
is suppressed. This behaviour is significantly different from processes with massive two- 
particle final states, where the phase space near threshold scales oc (3 and where soft-gluon 
radiation provides a source of large NLO corrections. 

At high energies the NLO partonic cross sections in the gluon-gluon and (anti)quark- 
gluon channels approach non-zero limits (plateaus) asymptotically, resulting in a ratio 
f<w /fgq 2Ca/Cf = 9/2. This has to be contrasted with the effective scaling behaviour 
oc 1/s of the LO cross section in the invariant energy range of the partons far above 
threshold, v/§ /i 0 (but at the same time hi y/§ hi m t , so that f-channel top-exchange 

diagrams are not yet logarithmically enhanced). The enhancement of the NLO diagrams 
is caused by the exchange of soft, nearly on-shell, space-like gluons, associated with the 
radiation of a hard parton from the initial state (see Fig. ||j). This gives rise to a radiative- 
return phenomenon, causing the gg-fusion subprocess induced by the virtual gluon to 
effectively take place at much lower energies, i.e. \/seff = 0(/i 0 ). I 11 analogy with the 
Drcll-Yan process, the cross section for the qq channel vanishes asymptotically in leading 
and next-to-leading order. Therefore we expect the hadronic cross sections at the Tevatron 
to be affected substantially less by these high-energy plateaus than the corresponding cross 
sections at the LHC. 

The average CM energy (V§) for the parton processes at the Tevatron and LHC 
colliders is approximately 650 GeV and 850 GeV, respectively, corresponding to 77 « 1 — 2 . 
From Fig. |5] we can conclude that in this range of 77 the NLO corrections at the central 
scale po are moderately negative for the quark-antiquark channel, which dominates the 
cross section at the Tevatron, and moderately positive for the gluon-gluon channel, which 
is the prominent production mechanism at the LHC. The (anti)quark-gluon channel is 
strongly suppressed at both the Tevatron and the LHC. 


3.2 Hadronic cross sections 

The numerical analyses of the hadronic cross sections have been performed for the 
Fermilab Tevatron pp collider with a centre-of-mass energy of y/s = 2 TeV, and for the 
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Figure 5: The scaling functions f a b for the partonic processes ab —> ttH as a function of 
r) = s/ 4yUg — 1. The notation follows Eq. (3T). The inserts show the ratio of the analytic 
threshold approximation and the complete numerical calculation. The Higgs-boson mass 
has been set to M H = 120 GeV. 
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9,9,9 

Figure 6 : Generic diagram leading to the high-energy plateau in the scaling f un ctions 

4/-4/ and 4M4 

CERN LHC with a pp centre-of-mass energy of y/s = 14 TeV. The hadronic cross sec¬ 
tions are obtained by convoluting the partonic cross sections with the parton distribution 
functions of the initial-state hadrons as specified in Eq. (|1 . 2|) . 


3.2.1 Total cross section and scale dependence 

The total cross section at the Tevatron is displayed in Figs. [7| and |[ Representative 
results are listed in Table [lj. For a Higgs mass between 100 and 150 GeV, the cross section 
varies between about 10 fb and 1 fb, if the central value p —> po = ( 2 mt + Mh )/2 is chosen 
for the renormalization and factorization scales. In NLO the theoretical prediction is 
remarkably stable with very little variation for p between ~ /i 0 /3 and ~ 3po, in contrast 
with the Born approximation, for which the production cross section changes by more than 
a factor of 2 within the same interval. It has been checked that no accidental compensation 
of scale dependences is introduced by identifying renormalization and factorization scales. 
Although the cross section at the Tevatron is strongly dominated by the qq annihilation 
channel for scales p ~ po, the proper study of the scale dependence requires the inclusion 
of the gg, qg, and gq channels. If p is chosen too low, large logarithmic corrections 
spoil the convergence of perturbation theory, and the NLO cross section would even turn 
negative for p < po/5. 

As apparent from Fig. | 8 |, the Ii factor, K = ctnlo/o'lo with all quantities calculated 
consistently in lowest and next-to-leading order, varies from ~ 0.8 at the central scale 
p = p 0 to ~ 1.0 at the threshold scale p = 2 p 0 . The K factor is nearly independent of 
Mh in the relevant Higgs-mass range (see Table p]). As argued in Ref. [0, the K factor 
below unity is in qualitative agreement with results obtained in the fragmentation picture 
proposed in Ref. [ [48| . 

The cross section at the LHC is dominated by the gluon-gluon process, which receives 
positive NLO corrections in the vicinity of the central renormalization/factorization scale 
p 0 . At p 0 we obtain K ~ 1.2, increasing to ~ 1.4 at the threshold value p = 2 p 0 , 
see Figs. |9] and [10|. Again these values are nearly independent of Mh in the relevant 
Higgs-mass range and compatible at the qualitative level with estimates obtained in the 

m 


fragmentation picture 


The NLO corrections significantly reduce the renormalization 


and factorization scale dependence and stabilize the theoretical predictions for the cross 
section at the LHC. 

As argued in Section pT7T| . the impact of the Coulomb singularity and the threshold 
logarithms from soft-gluon radiation is suppressed because of the (5 4 behaviour of the 
massive three-particle phase space near threshold. Nevertheless the threshold region gives 
us insight into the different colour factors that determine the NLO corrections at the 
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M h [GeV] 

MRST(LO/NLO) 

CTEQ6 


OLo [fb] 

0 NLO [fb] 

I< 

qq : gg (LO) 

onlo [fb] 


120 

5.846(2) 

4.857(8) 

0.83 

0.99 

0.01 

4.871(8) 

Tevatron 

140 

3.551(1) 

2.925(4) 

0.82 

0.99 

0.01 

2.931(5) 

160 

2.205(1) 

1.806(2) 

0.82 

0.99 

0.01 

1.808(3) 


180 

1.393(1) 

1.132(1) 

0.81 

0.996 

0.004 

1.133(2) 


120 

577.3(4) 

701.5(18) 

1.22 

0.30 

0.70 

665.0(19) 

LHC 

140 

373.4(3) 

452.3(12) 

1.21 

0.33 

0.67 

431.8(13) 

160 

251.6(2) 

305.6(8) 

1.21 

0.35 

0.65 

290.6(8) 


180 

176.0(1) 

214.0(6) 

1.22 

0.36 

0.64 

203.4(5) 


Table 1: Total cross sections and K factors for pp —> ttH at the Tevatron and for 
pp —> ttH at the LHC. The renormalization and factorization scales have been set to 
p = /i 0 = (2 m t + M h )/ 2. The MRST [ [47] parton densities have been adopted as default, 
whereas in the rightmost column a comparison is performed with the CTEQ6 [ |49| parton 
densities. The relative weight of qq and gg initial states at LO is given in the second-to-last 
column. 


Tevatron and the LHC. Just like the partonic high-energy region, where the fragmentation 
picture applies and the high-energy plateaus dominate the NLO corrections, also the low- 
energy threshold region hints at K factors at the LHC that are substantially larger than 
the ones at the Tevatron (see Section [2.4| ). 

We have also studied the uncertainty in the cross-section prediction due to the error in 
the parametrization of the parton densities. To this end we have compared the NLO cross 
section evaluated using the default MRST [[47] parametrization with the cross section eval¬ 


uated using the CTEQ6 [ fE]] parametrization. The results are collected in the rightmost 
column of Table [3]. Both MRST and CTEQ include a systematic treatment of parton- 
distribution-function uncertainties and provide the means for a quantitative estimate of 
the corresponding uncertainties in the cross sections. Using the parton-distribution-error 
packages [ |47|, |49|], we find that the uncertainty due to the parametrization of the parton 
densities in NLO is less than ps 5% at the Tevatron, where the quark-antiquark initial 
states dominate, and less than ps 10% at the LHC, where gluon-gluon initial states are 
more prominent. 


3.2.2 Differential distributions 


In Figs. p~THl5l the normalized transverse-momentum and rapidity distributions of the 
Higgs boson and the top quarks [t and t not discriminated] are shown for pp —> ttH at 
the Tevatron in LO and NLO approximation. For the transverse-momentum distributions 
we choose the transverse mass as the default renormalization and factorization scale, i.e. 
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p 2 = Pth + an d H — Prt/i + m ‘t f° r Higgs-boson and top/antitop distributions, 
respectively, which provides a more natural choice for large transverse momenta.]] 

The scale dependence of the distributions in rapidity and transverse momentum of the 
Higgs boson and top quarks is greatly reduced at NLO. The improvement is exemplified in 
Fig- O which compares the transverse-momentum distribution of the Higgs boson, once 
calculated with the renormalization and factorization scales set to the transverse mass, 
and once with both scales fixed to /i 0 . We observe differences of up to 30% at LO, while 
the shape of the Pt,h distribution at NLO is practically identical for the two choices of 
scales. 

The ratio of the (normalized) NLO and LO transverse-momentum distributions of the 
Higgs boson, displayed in the insert in Fig. [T2|. reveals that the simple rescaling of the 
LO Pt,h distribution by a constant K factor reproduces the NLO distribution only within 
±10% in the relevant transverse-momentum range, if the transverse mass is chosen as the 
scale. The increase of the ratio d<7NLo/duLO with increasing p T H can be understood in 
terms of the scale dependences of d<7NLO an d dcr LO : with p 2 = p^ H ± Mjj rising with 
increasing Pt,Hi doro decreases with increasing p faster than d<7NLO (see Fig. |8|). We note 
that choosing a fixed scale for the LO Pt,h distribution provides a better description of 
the NLO shape. 

For the transverse-momentum distribution of the top quarks, the increase of the ratio 
dc nlo/ dcr lo with increasing p T H due to the scale dependence is balanced by gluon radia¬ 
tion off top quarks at NLO, which reduces the energies and momenta of the top particles 
and thereby depletes the region of large transverse momenta. As a net result, the shape 
of the transverse-momentum distribution of the top quarks is barely affected by the NLO 
corrections in the relevant range pr,t/t < 300 GeV, if the transverse top-quark mass is 
adopted for the renormalization and factorization scales (see Fig. 14). 

The reduction in energies and momenta of the massive final-state particles due to 
gluon radiation at NLO is also expected to enhance the central rapidity region. However, 
as evident from Figs. [13] and [15], the effect on the (normalized) rapidity distributions of 
the massive final-state particles is small and becomes noticeable at the edges of phase 
space only. 

Cum grano satis , the transverse-momentum and rapidity distributions of the Higgs bo¬ 
son and the top quarks for pp —> ttH at the LHC, Figs. show the same qualitative 

behaviour as observed for the Tevatron. Adopting the Higgs transverse mass as the scale, 
the difference between the NLO and LO Higgs-boson transverse-momentum distributions 
at the LHC is naturally more pronounced, with deviations of up to 30% at large Pt,h- 
As for the Tevatron, choosing a fixed scale for the LO Pt,h distribution provides a better 
description of the NLO shape. 


4 Summary 

In this report we have presented the theoretical analysis of the Standard-Model pro¬ 
cesses pp/pp —> ttH ± X in next-to-leading order QCD at the Tevatron and the LHC. We 
have focussed on two aspects. 

5 The integrated cross sections a to t = f d px dcr/dpr reproduce the total cross sections with renormal¬ 
ization and factorization scales fixed to /x 0 within better than 5% at NLO. 
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(i) The technical basis for the prediction of the cross sections has been elaborated 
in detail. Particular emphasis has been put on evaluating pentagon diagrams, which 
involves the isolation of the singularities, the reduction to lower n-point functions, and 
the stable integration over the phase space. As a second important technical element we 
have applied the dipole subtraction formalism to treat massive final-state particles. In 
fact, the processes reported on here, represent the first complex examples that have been 
treated in this formalism in the case of massive quarks. 

(ii) In expanding on results presented in an earlier letter, Ref. [0. we have analyzed 
not only the total cross sections for Higgs bremsstrahlung off top quarks at the Tevatron 
and the LHC, but also final-state rapidity and transverse-momentum distributions for the 
Higgs boson and the top quarks. Moreover, we have studied the uncertainty in the next- 
to-leading order cross section prediction due to different parametrizations of the parton 
densities. 

When including higher-order QCD corrections, the K factor at the Tevatron has a 
value slightly below unity. It varies between ~ 0.8 and ~ 1.0 for renormalization and 
factorization scales between /i = /i 0 and 2/i 0 , almost independent of the Higgs-boson mass 
in the intermediate range. Here 2/i 0 = 2 m t + Mu denotes the threshold CM energy of 
the partonic subprocesses. Similarly, the K factor varies between ~ 1.2 and ~ 1.4 for 
the same scales at the LHC. Most importantly, the NLO predictions for the total cross 
sections and for the distributions in rapidity and transverse momentum of the Higgs boson 
and top quarks are stable when the renormalization and factorization scales are varied, 
in contrast to the Born approximation. The improved cross sections can therefore serve 
as a solid base for experimental analyses at the Tevatron and the LHC. 
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Figure 11: Normalized transverse-momentum distribution [GeV~ l ] of the Higgs boson for 
pp —> ttH+X at the Tevatron in LO (upper figure) and NLO (lower figure) approximation, 
for two choices of the renormalization and factorization scales: fi 2 = /j,q = (■ rrit + M^/2) 2 
and pi 2 = p^ H + Mg. 
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Figure 12: Normalized transverse-momentum distribution [GeV~ l ] of the Higgs boson for 
pp —» ttH + X at the Tevatron in LO and NLO approximation, with the renormalization 
and factorization scales set to p 2 = Pth + Mh- 



y H 


Figure 13: Normalized rapidity distribution of the Higgs boson for pp —> ttH + X at 
the Tevatron in LO and NLO approximation, with the renormalization and factorization 
scales set to po = (2m* + Mh)/ 2 . 
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Figure 14: Normalized transverse-momentum distribution [GeV~ l ] of the top and antitop 
quarks, (da/dpr,t + dcr/dpT,t)/2, for pp —> ttH + X at the Tevatron in LO and NLO 
approximation, with the renormalization and factorization scales set to p 2 = p^t/t + m ‘t ■ 



y t /i 


Figure 15: Normalized rapidity distribution of the top and antitop quarks, ( da/dy t + 
dcr/dyt)/‘2, for pp —> ttH + X at the Tevatron in LO and NLO approximation, with the 
renormalization and factorization scales set to po = (2 m t + Mjj)/2. 
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Figure 16: Normalized transverse-momentum distribution [GeV~ l \ of the Higgs boson 
for pp —> ttH + X at the LHC in LO and NLO approximation, with the renormalization 
and factorization scales set to p 2 = p^ H + Mjj. 



Figure 17: Normalized rapidity distribution of the Higgs boson for pp —> ttH + X at the 
LHC in LO and NLO approximation, with the renormalization and factorization scales 
set to no = (2 m t + Mh)/ 2. 
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Figure 18: Normalized transverse-momentum distribution [GeV~ l \ of the top and anti¬ 
top quarks, (da/dpT,t + da/dpr,?)/2, for pp —> ttH + X at the LHC in LO and NLO 
approximation, with the renormalization and factorization scales set to p 2 = p^ t /i + m 2 . 



y t /i 


Figure 19: Normalized rapidity distribution of the top and antitop quarks, (da/dy t + 
dcr/dyt)/‘2, for pp —> ttH + X at the LHC in LO and NLO approximation, with the 
renormalization and factorization scales set to po = (2 m t + 
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Appendix 

A The IR-divergent scalar 3- and 4-point integrals 


In this first appendix all relevant IR-divergent scalar 3- and 4-point integrals are 
listed. They have been derived in two independent ways, once with the help of Feynman 
parametrization, once by means of cutting techniques and dispersion integrals. In the 
following the momenta p t are defined as in Section 2.1 , i.e. they are on shell: p\ — p\ — 0, 
P 3 — Pa — m h P 5 = Throughout this appendix we make extensive use of the 

quantities 

' ” Pr ~ 1 


Pr = \ 1 


4m? 


x r = 


(A.l) 


r /3 r +1 

where r represents a generic invariant. Moreover, we adopt the convention that invariants 
with a bar receive an infinitesimally small, positive imaginary part: f = r + iO. 

There are four different IR-divergent 3-point functions that become relevant. In R ^ 4 
dimensions they read 


Co D) (Pi,~P 2 , 0,0,0) = ~ 


Co D \pi:Pl -*>3,0,0,771*) = 


a 2 K (r) - A i R (/4 ln ( + \ ln2 


nrit 
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4|) + u44| 
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■S 45 — ^23 


In 2 ( 1 


AlVjlnf 

S45 


V S45 — mi ) \rrit ) 


mi 


2 ,-ln 2 (i-^| 

mi J \ mi 


Cq D \p 3 , —Pi, 0 , m t , mf) = 


S. 34 R 


*34 


\ (p)H'Xs 34 ) -2U 2 (-x S34 ) 

1 7T 2 ' 

2ln(ats 34 ) ln(l + x, 34 ) + - ln 2 ^) - — 


(A.2) 


We also present these 3-point functions with mass regularization, where all internal 
massless particles are endowed with an infinitesimally small mass A: 


c, 


(mass,D=4) 

0 


(Pi,-P2, A, A, A) = 


2 s 


hr 


A 2 


G (mass,.D=4) 
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1 7T 2 
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(A.3) 

Taking into account that the decomposition of the 5-point integrals leads to additional 
4-point integrals, six different IR-divergent 4-point integrals are needed: 
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D { 0 D) {p 3,P3 - Pi,-pA,Q,mt,rnt,m t ) = 
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D ( 0 D) (p 3,P3 +P5, -P4,0,m t ,m t ,m t ) = 
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- 2 ln(z S34 ) ln^l - y§) - In 2 (xg) - ln 2 ( x m*) + y 


mi 

The following specific combinations of logarithms and dilogarithms were used, 

£i 2 (a;, y ) = Li 2 (l — xy ) + ln(l — xy ) ln(xy) — ln(x) — ln(y) , 

£i 2 (x, y, z) = Li 2 (l — xyz) + ln(l — a;?/*) ln(xyz) — ln(x) — ln(y) — ln(z) . (A.5) 

The expressions given for the C$ and D 0 functions are valid for arbitrary real values of 
the variables s, *jp and A/#, independent of their kinematical range. 

The method described in Section p.2.3| \(ii\ can also be used to calculate the singu¬ 
larities of the 4-point functions Dq in terms of 3-point functions Cq. The corresponding 
relations between Zt 0 , s ing and the 3-point subintegrals are 

D 0 (pi,Pi P3i P 2 j 0; 0, m t: 0) sing = -—-—j C 0 (pi,-P2,0,0, 0) 
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D 0 (P3iP3 + P5, -P4, o, m t , mt, m t ) s[ng = 


s 35 - m 


Cq(p3, —p 4 , 0, m t , m t ). 


(A.6) 


Obviously the explicit results for Dq D ' 1 and Cq 0 ' 1 given above are compatible with these 
relations. Moreover, the functions -Doling can be used to translate the D-dimensional Dq 
functions of Eq. ( |A.4|) into their mass-regularized counterparts D ( mass > D - 4 ), 
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(mass, D =4) 
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= D 


(mass,.D=4) 

0,sing 

(D) 
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n_ n( D ) 

u 0 -^O.sing 


where the C 0 functions occurring in D [ (j ' J J ng and D^ng ’° 4> are given in Eqs. ( |A.2| ) and 
(O), respectively, or result from those by mere substitutions. The dots represent in¬ 
finitesimally small terms proportional to the regulators (D — 4) or A. 

A simple example: 
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(A.8) 


which agrees with the result in the literature [ |29|. So, by subtracting and adding well- 
defined singular integrals we are able to switch between arbitrary regularization schemes. 
For instance, to switch to an off-shell regularization scheme, with some of the external 
particles put slightly off-shell, it is sufficient to know the relevant 3-point integrals in that 
scheme. 


B Colour operators and correlations 

In order to make the results in this report as compact as possible, colour indices have 
been suppressed in the expressions by making use of colour operators. The definitions of 
these colour operators are given in this appendix (see also Refs. [ 


24], H). 


Consider the matrix element for a process involving n coloured external particles, 
M ai - an , where oq,..., a n are the colour indices of the coloured particles. These colour 
indices range from 1 to 8 for gluons [adjoint representation] and from 1 to 3 for (anti)quarks 
[fundamental representation]. Emission of a gluon with colour index c from the external 
particle j is represented by a specific colour operator (charge) T, acting on the matrix 
element: 


(M 0 T,)«* = £ 4,,*, ■ • • (7}%* ■ • • S. nK M bl = £ (I}%* M 


cl\ ... bj ... dfi 


bi, , b n 


(B.l) 

with {Tf) ab = — if ca b for incoming/outgoing gluons, {Tf) ab = \ c ab /2 for outgoing quarks 
or incoming antiquarks, and {Tf) ab = — \ ba j2 for outgoing antiquarks or incoming quarks. 
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Owing to SU(3) gauge invariance, the complete matrix element .Ad is a colour singlet. 
This leads to the colour-conservation identity 

n 

®T i ) c;ai ■■■“" = 0. (B.2) 

3 = 1 


In the NLO calculations the colour operators enter as products (T,- • T*,) = (T*, • T ? ), 
dehned according to 

[M ® (T, . T,)]« •••- = E E (r/)a jtj (T t % ih M« ■■■"* ■■■“■ if j + k, 
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(B.3) 


Here Cj is the Casimir operator in the representation of particle j, i.e. = C\ = 3 
for gluons [adjoint representation] and T) = CV = 4/3 for (anti)quarks [fundamental 
representation]. At the cross-section level this gives rise to colour correlations of the form 


E lxil 2 ®(T r T t )= e 

colours ai,...,a n 




E E ( r /)«A Cr t % A M«- 1 ‘ 

bj 5 bjc c 

(B. 4 ) 

for j k. Needless to say that the left-hand side of Eq. ( |B.4| ) looks much more appealing 
from the point of view of compact expressions. 

Finally, we give the explicit results for the colour correlations in terms of the colour 
operators C qq and Cjjf defined in Eqs. (|2.6| ) and (|2.14|) , respectively. Since these C oper¬ 
ators represent a basis for the relevant colour structures of the qq and gg channels, the 
correlations C ® (T; • T*,) can again be expanded in terms of C’s: 


C«®(T,.TJ = $C«+..., C® ® (T, ■ T») = E 


(B.5) 


n =1 


where the dots stand for (irrelevant) operators orthogonal to C qq . The coefficients £ follow 
from simple colour algebra. For the qq channel they read 


1 7 l 

cqq — cqq — cqq — cqq — cqq — cqq — 

S12 — S34 — g> S13 — S24 ~ gj S14 ~ S23 — g • 

For gg fusion it is convenient to give them in matrix form, 

3 3 

£ gg,mn\ ’ ■ / ~ a a 


«??'”“) = diag (- 3 , - 5 ) , (&»•”“) = diag(-|, i,± 
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(B.7) 
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These £ coefficients have also been used in Section [2.2.5| in rewriting Eqs. (|2.45|) and (|2.46| ) 
in terms of Eq. ( |2.48| ). 
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C The Altarelli Parisi splitting functions at NLO 


The so-called Altarelli-Parisi splitting functions P ab (x) [ [5(J are a measure of the 
probability of finding the parton b with fraction x of longitudinal momentum inside the 
parent parton a|] At NLO these splitting functions can be written in terms of a regularized 
splitting function and an IR-sensitive part according to 


P a \x) = P T %(x) + 5' 


ab 


2 T 2 ' 1 

a ' 1-x 


+ 7a 5(1 -x) 


(C.l) 


The anomalous dimensions y a are given by 


7 9 = 7 9 -=2 C f = 2 


11 ^ 2 33 — 2N f 

79 ~ ~6° A 3 Tr f ~ 6 : 


(C.2) 


using the customary normalization Tr = 1/2. The ‘+’ distribution (..,) + is defined in 
the usual way: 


dx 


f(x) ,g(x)= / d xf(x) [g{x)-g{l)\ . 


(C.3) 


The regularized splitting functions are given by 


P&iz) 

PrlW 
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P r %(x) 
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p?U x ) 

PZ(x) 
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o, 


a 


C F (1 + x), 

1 + (1 — x) 2 
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T, 
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x 2 + (1 — x)' 


2 C\ 


T — x 
. x 


1 + x(l — x) 


(C.4) 


If the Altarelli-Parisi splitting functions were calculated in D — 4 — 2e dimensions 
by means of dimensional regularization, a term — eP' ab {;x) would have to be added to 
Eq. flC.lp . Such a term enters the final NLO expressions through the interference with 
the 1/e collinear poles. The associated splitting functions P' ab (x ) amount to 

Pqq( X ) = P'qq{x) = P' gg (x) = 0, 

Pqq(x) = PP(X) = Pqg(l~X ) = P' gg {l ~ x) = C'f(I- x), 

P n {x) = P g9 (?) = 2T r x(1-x). (C.5) 


Among other things, these splitting functions take into account the difference between the 
number of spin degrees of freedom of the quarks (2) and gluons (D — 2) in D dimensions. 


6 We have adopted the notation of Ref. [ , which may differ from the conventions elsewhere in the 

literature by the interchange of the parton labels a and b. 
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